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THEORETICAL SYMMETRIC SPAN LOADING DUE TO FLAP DEFLECTION FOR WINGS OF 

ARBITRARY PLAN FORM AT SUBSONIC SPEEDS 1 

By John DeYoung 


SUMMARY 

*4 simplified lifting-surface theory is applied to ike problem of 
evaluating span loading due to flap deflection ■ for arbitrary wing 
plan forms. With the resulting procedure, the effects of flap 
deflection on the span loading and associated aerodynamic 
characteristics can be easily computed for any wing which is 
symmetrical about the root chord and which has a straight quar- 
ter-chord line over the wing semispan. The effects of compressi- 
bility and spanwise variation of section lift-curve slope are 
taken into account by the procedure. 

For the case of straight-tapered wings, the load distribution 
and lift due to flap deflection are presented in chart form for a 
broad range of wing plan forms. Experimental values of flap 
effectiveness are compared with the values obtained from these 
charts for several wings. The comparison shows good agree- 
ment between theory and experiment. 

The method presented can also be used to calculate the down- 
wash in the vertical center of the wake of a wing which has 
arbitrary spanwise loading. 

INTRODUCTION 

Reference 1 presents a procedure which is based on simpli- 
fied lifting-surface theory and which permits the rapid 
determination of the symmetric spanwise distribution of lift 
on wings having essentially arbitrary plan form. With 
reference 1, the symmetric spanwise loading can be found for 
any symmetric angle-of-attaek distribution such as that due 
to spanwise continuous twist or camber. The method satis- 
fies the boundary condition at only a finite number of span- 
wise points; thus, to insure accuracy for the loading obtained 
by this method, the angle-of-attack distribution must be 
continuous and must not vary excessively. If the angle-of- 
attack distribution is discontinuous, such as that due to 
partial wing-span flaps, or varies excessively, the boundary 
condition should be satisfied at sufficient points to insure a 
correct answer. This normally entails an excessive amount 
of labor, since the computations required increase approxi- 
mately as the square of the number of points used. 

An alternative procedure is to use a method which satisfies 
the boundary condition at only a few points where this 
boundary condition is adjusted so that at these points the 
calculated loading is identical with that which would be 


‘ Supersedes NACA TN 2278, "Theoretical Symmetric Span Loading Due to Flap Deflec- 
tion for Wings of Arbitrary Plan Form at Subsonic Speeds” by John DeYoung, 1951. 


predicted by a method satisfying the boundary condition 
at all points. This specific boundary condition for the case 
of antisymmetrically deflected ailerons is developed in refer- 
ence 2, which is based on the zero-aspect-ratio theory of 
reference 3. A similar development can determine the 
specific boundary condition for the case of symmetrically 
deflected flaps. 

It is the purpose of the present report to determine the 
specific boundary condition and to provide simple methods 
for computing the aerodynamic characteristics for wings with 
symmetrically deflected flaps. 
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NOTATION 


aspect ratio 



coefficient depending on wing geometry and 
indicating the influence of symmetric load- 
ing at span station n on the downwash 
angle at span station v 

span of the wing measured perpendicular to 
the plane of symmetry, feet 
wing chord, 2 feet 


(f) 


mean wing chord 2 
local lift coefficient ( 
lift coefficient ( 


> feet 


/local lift 


V 


qc 

^ total lift\ 


qs 


rate of change of lift coefficient with flap 
deflection, per radian 

rate of change of lift coefficient with wing 
angle of attack, per radian 


spanwise loading coefficient for unit lift 


coefficient 



scale factor 

factors of loading interpolation function 
spanwise loading coefficient or dimensionless 


circulation 



spanwise loading coefficient per radian 


flap deflection 



of 


• Measured parallel to the plane of symmetry. 
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wing-chord distribution parameter (at span 
station v) which includes the effects of 
compressibility and section lift-curve slope 



integration factors for spanwise loading due 
to flaps 
Mach number 

arbitrary number of span stations 
free-stream dynamic pressure, pounds per 
square foot 
wing area, square feet 
free-stream velocity, feet per second 
induced velocity, normal to the lifting surface, 
positive for downwash, feet per second 
longitudinal coordinate measured from the 
lateral plane through the quarter chord of 
the wing-root chord , 3 feet 
lateral coordinate measured from the wing 
root perpendicular to the plane of sym- 
metry, feet 

wing angle of attack , 3 radians 
section angle of attack at span station v , s 
radians 


lift-effectiveness parameter 3 



compressibility parameter Vi 1 — M*\ 
circulation, feet squared per second 
angle of deflection of flap , 3 radians 
angle of deflection of flap measured perpen- 
dicular to the hinge line, radians 


dimensionless lateral coordinate 



dimensionless flap span on one wing panel, 
measured perpendicular to the plane of 
symmetry, from the wing root outboard 
for inboard flaps, and from the wing tip in- 
board for outboard flaps 

trigonometric spanwise coordinate defining 
the spanwise end of the flap, radians 
ratio of section lift-curve slope at span sta- 


tion v to — , both at the same Mach 


number 4 


sweep angle of the wing quarter-chord line, 
positive for sweepback, degrees 
sweep parameter which includes the effects of 

compressibility (tan -1 degr< 


;rees 


/ tip chord \ 

\root chord/ 

kk) 

trigonometric spanwise coordinate (cos 1 77) , 
radians 


wing taper ratio ( - 
dimensionless longitudinal coordinate 


* See footnote 2, page 1, 
4 See*footnote 2, page 1, 
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The simplified lifting-surface theory used herein replaces a 
lifting surface by a lifting vortex located at the wing quarter- 
chord line. The boundary condition specified for determin- 
ing the vortex strength distribution is that, along the three- 
quarter-chord line of the wing, there shall be no flow through 
the lifting surface. 

The mathematical method of solution and discussion of the 
limitations of the theory are given in detail in NACA Rep. 
921, 1948 (reference 1 ). The same mathematical procedures 
and assumptions are used in the present report. 

The mathematical solution for the determination of span- 
wise loading consists of satisfying the boundary condition 
that along the three-quarter-chord line of the wing there shall 
be no flow through the lifting surface. If the boundary 
condition (i. e., the spanwise twist along the three-quarter- 
chord line) can be approximated by a few terms of a mathe- 
matical series, then the “no-flow” condition need be satisfied 
exactly at only a few spanwise stations and the calculated 
loading can be expressed by a few terms of a mathematical 
series. Thus it has been found that satisfying the boundary 
conditions at four spanwise points enables accurate compu- 
tation of the loading due to symmetric twist (reference 1 ), 
whereas only three spanwise points are required for accurate 
computation of the loading due to antisymmetric twist 
(reference 2 ). The distribution of spanwise twists in these 
two cases is such that at. most a four and a three term trig- 
onometric series, respectively, are needed to duplicate them. 

For the case of irregular and discontinuous distribution of 
twist such as that required to simulate the span loading due 
to deflection of a partial-span flap, it becomes necessary to 
satisfy the boundary condition at a large number of points, 
particularly near the point of discontinuity, in order that 
sufficient accuracy will be insured. Since the computation 
labor increases approximately as the square of the number of 
points considered, the foregoing method becomes less prac- 
tical. An alternate approach is to alter the boundary con- 
dition for the method involving few points at which the 
boundary condition is satisfied, so that at these points the 
loading is identical with that obtained by the method for 
which the boundary condition is satisfied at all points. 
This process is practical because the spanwise loading of 


integer pertaining to span station associated 
with the loading (v= cos 1l ' n ' 


to 4 - 1 , 

integer pertaining to span station associated 
with the downwash ( w=eos 

V to + 1, 

half-integer pertaining to span station 

( ‘kir 

I 77 = cos 


m 4-1, 
pertaining to flaps 
wing tip, 
wing root 
average or mean 


denoting full wing-chord flaps 
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wings becomes independent of plan form as aspect ratio 
approaches zero (reference 3). Thus the solutions for the 
zero-aspect-ratio wing which satisfy the boundary condition 
at all points can be used to obtain the proper adjustment of 
the boundary condition for a method that satisfies the bound- 
ary condition at only a few points. With the altered bound- 
ary condition, the method satisfying conditions at only a 
few points will give accurate results for plan forms of very 
low aspect ratio, and should give good accuracy for wings 
of high aspect ratio since the alteration of the boundary 
condition can be kept small by proper choice of flap spans, 
and also since loading distribution is affected less by plan 
form than by the boundary condition. The mathematical 
procedure for determining the altered boundary conditions 
is developed in appendix A. 

Although an altered boundary condition can be found for 
an arbitrary position of the point of discontinuity, the 
alteration can be made a minimum and of small order by 
locating the discontinuity at the mean position between 
two points at wliich the boundary condition is satisfied. 
With the point of discontinuity defined, the flap span be- 
comes defined. The results for arbitrary flap spans are 
found by cross plots of results for the defined spans. 

When the flap spans and corresponding twists (altered 
boundary conditions) have been determined, the method is 
applied in the same manner as that of reference 1. In this 
way, (m-rl)/2 linear equations in terms of symmetric- 
loading distribution are obtained which satisfy the wing 
angle-of-attack conditions 5 at the three-quarter-chord line 
at m stations, where m is an arbitrary odd integer. The 
application in the present report is with m— 7. The linear 
equations are represented by the summations 

«,=2 7=1,2, 3, 4 (1) 


where the integer subscripts v and n pertain to span stations 
?? = cos (vt/S) and cos (titt/S), respectively. To obtain the 
loading coefficients G n from the simultaneous solution of 
equation (1), it remains only to evaluate the coefficients 
fl,* (which depend on wing geometry, compressibility, and 
section lift-curve slope) and the spanwise variation of the 
symmetric angle of attack a,. The reader should note that 
with a, constant (spanwise) equation (1) represents additional 
loading, and that with a, varying symmetrically and con- 
tinuously it represents basic loading (provided the net lift 
is zero). 

EVALUATION OF a„ 

As shown in reference 1, a,* can be conveniently presented 
as a function of two parameters, 6 namely, a sweep parameter 
defined as A^=tan -1 [(tan A) ,//3] and a wing-chord distribution 
parameter H, defined by 


//, = d. 



* As discussed in reference 1, a, is not limited to small values but am be as large an angle 
as desired, provided separation does not occur. 

* The effects of compressibility and section lift-curve slope are equivalent to a change in 
wing plan forms and can be taken into account by a proper adjustment of the a,*, values. 


where 

k, ratio of experimental section lift-curve slope 

at span station v to the theoretical value of 
2 T~ 

— > both at the same Mach number 
P 

c„ wing chord at span station v 

d, scale factor which has the following values: 

0.061 for v=l 
.234 for v=2 
.381 for r=3 
.320 for r=4 


Equation (2) can be written in an alternate form that gives 
H, in terms of wing geometry parameters that are more 
significant; thus 


where 


H,=d, 


d M ) 

V Kav / 


1 

(a'./ Vip) (£*/ Cep) 


( 3 ) 


A'ac 

K 

K<2B 


An 

C an 

£A 


ratio of average section lift-curve slope to 
27r//3 , both at the same Mach number 

spanwise distribution of section lift-curve 
slope for a given Mach number 

spanwise distribution of the wing chord 
aspect-ratio parameter 


The term 1 /'(k/0 (c,jc ae ) of equation (3) gives the effective 
aerodynamic taper of a wing. The distribution of *:/*„ may 
vary with Mach number, particularly at transonic speeds 
(e. g., due to spanwise variation of airfoil section). However, 
since this distribution contributes to taper effect, the primary 
change will be in the loading distribution rather than in gross 
load. 

With H y determined from equation (2) or (3), the values of 
o,», 16 in all, can be found in figure 1 where a,* is given as a 
function of H, for various values of As. 

For the case of straight-tapered wings with arbitrary sec- 
tion lift-curve-slope distribution for which the chord distri- 
bution is specified by taper ratio, evaluation of equation (3) 
is facilitated through the use of figure 2, which gives KyH,jfiA 
for each of the four span stations as a function of taper ratio. 


EVALUATION OF a,. 

The symmetric angle-of-attack distributions most com- 
monly encountered are those resulting from constant angle of 
attack, twist and camber, and the effective angle of attack 
due to flap deflection. The first two cases are the subjects 
of reference 1. The effective angle-of-attack distributions 
due to flap deflection are the subject of the present report. 

To minimize the computations involved, it is convenient 
to consider the cases of inboard and outboard flaps inde- 
pendently. The inboard flaps are measured from the wing 
midspan outboard, and the outboard flaps are measured 
from the wing tip inboard. Three flap spans for inboard 
and outboard flaps are considered in appendix A and the 
equivalent angle-of-attack distribution for unit flap deflection 
a y /S L for each case is listed in table A4. The values of a,/S L 
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given in table A4 are for (c f jc) = 1. The results can be easily 
modified for the case of (c,/c)<l as will be shown later. 
The case of full wing-span flaps is easily considered since 
this is the case of constant angle of attack which results in 
additional-type loading. 

APPLICATION OF METHOD 

For the cases of symmetric angle-of-attack distributions 
over a wing with flaps deflected, it is possible to form the set 
of simultaneous equations which are required for the solution 
of loading distribution for an arbitrary plan form. With the 
loading known, integration formulas can be given to deter- 
mine aerodynamic coefficients. 

The loading distribution G n , determined from the solutions 
of the simultaneous equations, is a function of a, n , which has 
been shown in a preceding section to be a function of wing 
geometry, compressibility, and section lift-curve slope. The 
aerodynamic coefficients are obtained by integrations of the 
laoding distribution and, therefore, will also be functions of 
wing geometry, compressibility, and section lift-curve-slope 
parameters. 

In the following sections, application of the method to 
obtain the general solution for arbitrary spanwise distribu- 
tion of wing chord is outlined and solutions are presented for 
the special case of wings having straight taper. 

GENERAL SOLUTION 

Simultaneous equations for the loading.— The a vn values 
to be substituted in the equations that follow are obtained 
from figure 1 and table 1 7 with values of H, given by equation 
(2) or (3). 

1. Deflected inboard flaps: The flap spans measured from 
the wing midspan outboard are given by i y. The set of 
simultaneous equations for symmetric spanwise loading due 
to deflection of any of the tliree following flap spans can be 
obtained from the appropriate set of the following equations: 


Case 

I 

II 

m 

Vf 

0.195 

0. 556 

0.831 

ai 

Si 

-0.017 

0. 032 

-0.041 

<32 

Si 

.029 

-.003 

1.021 

<33 

5l 

-.014 

.994 

.955 

<34 

Si 

.98S 

.976 

1. 010 


= Q2 1 ^l-l" C?2 T 0 23 ^3 + 02t 

>- 

= G31 C32 5*2 -f <*33 5*3 +<*34 Gi 


where 


Gn—GJhi 


(4) 


For the case of ij/=l ; a v /bi is unity over the span of the wing 
(same condition as solution for. additional loading). 

2. Deflected outboard flaps: The flap spans measured from 
the wing tip inboard are given by ?? s . The se t of simultaneous 
equations for symmetric loading due to deflection of any of 
the three following flap spans can be obtained from the ap- 
propriate set of the following equations : 


7 Values of a 9n beyond the scope of figure 1 and table I ean be obtained from equation (B2), 
which gives the linear asymptotes of the a vn function. 


Case 

IV 

V 

VI 



0.169 

0.444 

0.805 


<31 

Si 

1.041 

0.96S 

1.017 

** ou Sr -f-c \iGi -f-<7 13 Oi ■¥ d 

ccj 

Si 

-.021 

1.003 

.972 

* cz-21 G i -j" o 22 Gi -f tt23 Gi 4- a j \ G\ 

03 

Si 

.045 

.006 

1.014 

“ a 3i Gi -f-aa2<?2 -f d 33 Gii-an G\ 

Cti 

Si 

-.040 

.024 

.012 

= cn CJi-f-c 12 Oa -f-au G% -f a u(7* 


( 5 ) 


where 

For the case of ij^=l, afbi= 1. 

Loading distribution. — The spanwise loading distributions 
can be determined for various flap configurations as follows: 

1. Full wing-chord flaps \{c f jc) = 1]: The loading is zero at 
the wing tip and is determined at four other span stations of 
??=cos j-nir/8) =0.924, 0.707, 0.383, and 0 by solution of the 
appropriate set of simultaneous equations from the previous 
section. With equation (A9) and tables A6, A7, and A8, 
the loading can be interpolated at span stations ij=0.981, 
0.831, 0.556, and 0.195 for eacli of the flap spans considered. 
With these given points and the knowledge that the slope 
of the loading-distribution curve is theoretically infinite at 
the point of angle-of-attack discontinuity or flap spanwise 
end, the loading distribution can be faired. For flap spans 
other than those given, the loading distribution of flaps that 
are neither inboard nor outboard can be obtained from the 
concept that loadings are additive (discussed later). 

2. Constant fraction of wing-chord flaps [(cflc)= constant]: 

The spanwise loading distribution is equal to the product of 
the loading distribution due to full wing-chord flaps and the 
lift-effectiveness parameter dajdb. The factor dajdb is a 
function of the ratio of flap chord to wing chord c f jc. Values 
of dajdb are presented in figure 3, which is reproduced from 
figure 18 of reference 4. Although reference 4 from which 
figure 3j\vas taken limits the Mach number range to values 
less than. 0.2, this limitation is believed to be too restrictive 
since linearized theory indicates that dajdb is unaffected by 
compressibility for the two-dimensional wing. As indicated 
in reference 3, dajdb is strongly affected by aspect ratio if 
the aspect ratio (or / L4 ) becomes very small; however, com- 
parisons with experiment indicate the values of dajdb from 
figure 3 appear to be valid for /3A> 2. This limit to which 
dajdb can be used will be considered further in the section 
Discussion when comparisons with experimental results are 
made. _ ___ 

3. Arbitary spanwise distribution of flap chord [( c f jc) = 
variable \: The flap can be divided into several parts each 
having constant dajdb and the load distribution due to each 
part determined. The total distribution is then the sum of 
these individual load distributions. 

Lift coefficient. — The lift coefficients due to flap deflection 
for the same flap configurations previously discussed can be 
found as follows : 

1. Full wing-chord flaps [{c,jc) = 1]: The spanwise loading 
due to flap deflection is, in general, so complicated that 
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numerical integration based on four spamvise values cannot 
be performed accurately with conventional integration 
formulas. However, in appendix A a special integration 
formula (which needs but four spamvise values) is developed 
which applies to the different flap spans. Equation (A6) 
can be written as - 

— -=(— )(h l G l +k&+h 3 8 a +k i G i ) ( 6 ) 

where, for each of the cases of equations (4) and (5), the h n 
values are given by 



Inboard flaps 

Outboard flaps 

, Fun 
span 
flaps 

Case 

I 

n 

m 

IV ' 

V 

VI 

V 

0.195 

0.556 

0.831 

0.169 

0. 444 

0.S05 

L 000 

hi 

0.299 

0. 299 

a 300 

0.302 

0. 301 

0.301 

0.301 

h; 

.544 

. 544 

.555 

.561 

.556 

.556 

.555 

h 

.725 

.725 

.725 

.734 

.728 - 

.726 

.726 

h , 

.392 

.392 

.392 

.405 

.395 

.393 

.393 


2. Constant fraction of wing-chord flaps [(c//c) = constant]: 
The flap effectiveness is given by 


&C l s da A 

Kqs d5 \ Kq 3 J 


( 7 ) 


3. Arbitrary spanmse distribution of flap chord [(c//c) = 
variable): Flaps for which cfc varies spanwise on the wing 
can be considered as equivalent to a wing-twist distribution. 
The effective symmetric twist of the wing is given by 


da . 
a ’ = d5 5 


(8) 


where da/db is now a function of spamvise position. If 
equation (8) is a continuous distribution such, that it can be 
plotted by specifying its value at four semispan points, then 
the wing can be considered to be twisted and solutions found 
directly as for basic loading. 

When equation (8) is discontinuous, the angle of attack 
can be divided into spanwise steps of constant angle of 
attack and the total lift can be found by the summation of 
the lift due to each spanwise step. The lift of a spanwise 


step is obtained from a curve of lift coefficient 





as a 


function of the extent of flap span from the wing root out- 
board. Such a curve can be obtained from the solutions. of 
equation (6). As previously indicated for r\ s =l 



The lift coefficient, due to the twist given by equation (8) 
can be obtained by a method other than the step method 
from the integral given by 


dflM in ( 9 ) 

Jo dS dr; 

which can be integrated numerically by taking the graphical 
slopes of pCujjKat. 

C724S3 — 54 17 


Spanwise center of pressure, induced drag, and pitching 
moment. — The integration formulas of loading due to flaps 
required to obtain spanwise center of pressure and induced 
drag are not developed here. However, the integration 
formulas given in reference 1 to obtain these characteristics 
can be used with acceptable accuracy for the ease of loading 
due to flap deflection. 

No rigorous procedure for determining pitching moments 
can be developed from this method. As noted in reference 
1, however, a good approximation for a wing without flaps 
can be obtained if the wing sweep is large since moments due 
to shifts of spanwise load overshadow those due to shifts of 
chordwise load. Where the effect of flaps is to be considered, 
as in the present case, it appears unlikely that the shifts in 
chordwise loads can be safely neglected. Therefore an 
estimation of pitching moments due to flaps should not only 
include the moment due to spanwise load redistribution, 
which this method will give accurately from the longitudinal 
moment of the center of load, but also should include an 
estimation of the moment due to chordwise redistribution. 

Additional considerations. — Several items pertaining to 
the usage of the method should be considered. 

1. Additice nature of loading and spanmse angle of attack: 
The linear relation between angle of attack and loading 
distribution of equation (1) states that all loadings are 
additive if the respective angle-of-attack distributions are 
added. 

Thus, the aerodynamic coefficients that result from the 
integration of the loading distribution are additive if they 
depend linearly on the loading distribution. Hence, lift 
coefficient, rolling-moment coefficient, and pitching moment 
are additive; whereas induced drag, spanwise center of 
pressure, and aerodynamic center are not. 

The additive concept is very useful for the determination 
of loading due to flap deflection. The loading due to an 
arbitrary span flap having arbitrary position on the wing 
semispan is simply found by adding or subtracting known 
loadings due to flaps at other locations on the wing. For 
example, if the spanwise ends of the flaps are at the same 
span stations, then for full-chord flaps 


G inboard — ' G T .^ \ ffoutboard ") 

( ^ujlnbosrd ( ^ i ^)outboard 

^ ^outboard 


( 10 ) 


2. Differential flap angles: The effect of a differential 
between flap angles can be taken into account by considering 
C Ll of each wing panel as having one-half the value of the 
symmetric case. The total wing lift is then the sum of the 
products of C L f2 and the angle of deflection of each flap. 
Although the total wing lift can be found by this procedure, 
the loading distribution cannot be found exactly by the 
products of the symmetric G/S and the deflection of each flap. 
However, the loading distribution so obtained will be suf- 
ficiently accurate for most purposes since this procedure 
neglects only the small change due to the induced effects of 
the difference in flap deflection on each opposite wing panel. 

3. Flap angles measured perpendicular to the hinge line: 
The relationship between a flap angle measured in a plane 
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perpendicular to the flap hinge line and that measured in 
a plane parallel to the plane of symmetry is given by 


tan S 


tan 8 

COS Ay 


( 11 ) 


where Ay is the sweep angle of the flap hinge line, and 8 is the 
angle measured perpendicular to the hinge line. Equation 
(11) can be inserted into any of the equations involving 8 to 
put the equations in terms of 8. 

For constant fraction of wing-chord flaps on straight- 
tapered wings, Ay is given by 


tan Ay=tan A 


4[(3/4) — ( C/ /c)] / l— X \ 

-1 Vl+X ) 


( 12 ) 


SOLUTIONS FOR STRAIGHT-TAPERED WINGS 

Charts of aerodynamic characteristics for straight-tapered 
wings -can be presented in terms of parameters involving 
geometry, compressibility, and section lift-curve slope. 8 
These charts provide a ready means of obtaining data 
directly. 

Equation (4) has been used to evaluate the loading dis- 
tribution (?/$i = (1/Si) (cic/2b) for a range of wing parameters 
and flap spans, and the results are plotted in figure 4. The 
loading presented in figure 4 gives the values due to full 
wing-chord flaps. For other flap configurations, the pro- 
cedure given in the loading distribution section of the general 
solution can be used. 

With the loading distribution determined from equation 
(4), the application of equation (6) provides lift coefficients 
due to flap deflection in the form PC Li Jk„„ for several flap 
spans. In figure 5, these values are plotted against extent 
of flap span from the wing root outboard for a range of wing 
parameters. 

As presented, figure 5 gives directly the lift due to full 
wing-chord inboard flaps for flap spans measured from the 
plane of symmetry outboard. The lift due to outboard 
flaps for flap spans measured from the wing tip inboard can 
be obtained from figure 5 by use of the relations of equation 
(10). For full wing-chord flaps located arbitrarily on the 
wing semispan, the lift can he obtained from figure 5 as in- 
dicated in the following example sketch: 



' Throughout the figures, is the constant spanwise-seetfon lift-curve slope dr the average 

of a small variation. For large span-wise variations of « that follow the function given in equa- 
tion (B5) developed in appendix B, the parameters and A can be replaced by the 

0/1(1+ A) , , . 

parameters > , and - A, respectively. 

Or+*<) Kr 

For large span wise variations of k that do not follow' the curve of equation (B5), the simul- 
taneousequationsfor the general solution can besolved for arbitrary distributions of «, Values 
of Hw can be obtained conveniently from figure 2. 


With the full wing-chord values given above, the lift due 
to constant fraction of wing-chord and flaps of arbitrary 
spanwise chord distribution can be found through use of 
equations (7) and (9) with the dajdS values of figure 3. 

APPLICATION OF THE a, n COEFFICIENTS TO THE PREDICTION 
OF DOWNWASH FROM ARBITRARY WINGS 

In the subject report and in reference 1; it is shown how 
various forms of symmetric loading can be found quickly 
through use of simplified lifting-surface theory and values of 
the a vn coefficients which are presented in graphical form. 
It can also be shown that if the symmetric load distribution 
is known, from any computational method or from experi- 
ment,- then the same a, n coefficients can he used to find cor- 
responding values of the downwash in the center of the down- 
wash wake. It is the purpose of this section of the report 
to make available the method for computing downwash 
without attempting to explore the. possibilities or limitations 
of the method. In view of the importance of being able to 
predict downwash, it is considered that inclusion of this mate- 
rial herein is justifiable although not directly related to the 
main purpose of the report. It is important t‘o remember 
that the method can be used for any case of symmetric load- 
ing and is not restricted to that discussed in this report. 

As noted previously, the a vn coefficients provide a very 
simple means of finding the downwash in the center of the 
downwash wake of the wing due to any symmetric spanwise 
loading distribution for any sweep angle. The downwash 
can be found at specific points in the region between the 
wing tips from near the quarter-chord line to an infinite 
distance downstream. 

Equation (1) was derived (reference 1) to find the down- 
wash. at the effective tliree-quarter-chord line for which the 
distance downstream from the quarter-chord line is k(c/ 2). 
The general expression for the distance downstream from 
the quarter-chord line is 

x—jSy tan A? 

where x is the longitudinal distance measured in the hori- 
zontal plane from the wing-root quarter-chord point to an 
arbitrary downwash point. Substituting this general ex- 
pression. in equation (2) gives 

//,= 7V7- „ r — T 7- - (13) 

(w) — ij f tanA fl 

where 

4 6/2 

Equation (1) (for values of the parameters £//J, rj„ and tan 
A#) gives the downwash angle as 

(W) f = S a "n G„ (14) 

where the a tri coefficients are obtained from figure 1 with the 
H„ values given by equation (13). 

It should be noticed that, because the chordwisc loading 
is concentrated at the one-quarter-chord line, equation (13) 
is independent of aspect ratio and spanwise distribution of 
wing chord, while from equation (14) it can be seen that 
downwash depends on these wing geometric values only 
through the loading distribution G n . 
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Equation (14) 'will give the downwash angle due to arbi- 
trary symmetric' loading distribution for all values of 

P 

tan A;}. The downwash can be found for the specific span 
stations corresponding to v=l, 2, 3, 4, or ij,= cos -^-=0.924, 

O 

0.707, 0.383, and 0, respectively. If compressibility effects 
are included, the loading distribution must be that at the 
Mach, number considered. 

For convenience, equation (14) is rewritten such that 
downwash angle and loading distribution are in terms of 
unit lift coefficient. Then 


( w/V A = 1 

\pc z jK tt Ji tVr 9 (§A'\ 



(15) 


where c t cjC L c at is the span loading distribution for unit lift. 
Both a rn and c t c!C L c a , are affected by compressibility, being 
functions of the parameters pAjiCa, and A$. 

Equation (14) or (15) gives the downwash, in the vertical 
center of the wake, due to the integrated effect of the loading 
of a swept vortex and the trailing vortex sheet extending 
'backwards to infinity. Any rolling of the trailing vortex 
sheet into two trails is -neglected. Since either of the equa- 
tions are for arbitrary loading distribution, the loading 
coefficient can be set equal to a constant to duplicate uni- 
form loading and the downwash results obtained for the 
case of the trailing sheet rolled into two trails. 

With c t c/C L c a ,— 1 , equation (15) becomes 



and b r (as used in the past from line-theory concepts) is the 
span between the two trails given by the integrated loading 
divided by the loading at the middle of the wing, or 


fc ? = o 


DISCUSSION 

INFLUENCE OF WING GEOMETRY ON FLAP EFFECTIVENESS FOR 
STRAIGHT-TAPERED WINGS 

To show the effect of wing geometry and compressibility 
on lift due to inboard flap deflection, the data of figure 5 
have been cross-plotted and presented in figure 6. These 
effects are similar to those on wing lift-curve slope shown in 
reference 1. Figure 6 (a) compares the present theory with 
that of reference 3 and shows, as might be expected, that 

B 4. 

both theories give essentially the same results for 

Figure 6 (b) shows the influence of taper ratio on flap effec- 
tiveness for swept wings. Figure 6 (c) gives the relative 


effects of the wing parameters of sweep, aspect ratio, and 
taper ratio. 

Figure 7 shows comparative effects of inboard and outboard 
flaps on swept wings and shows the loss of flap effectiveness 
due to sweepback is less for inboard flaps. The relative 
effects indicated by figures 6 and 7 apply equally well for 
constant fraction of wing-chord flaps since the lift is only re- 
duced by a constant spamvise factor da/d8. 

LIMITATIONS OF THE dajdS CONCEPT AND COMPRESSIBILITY 
CONSIDERATIONS 

A minimum value of aspect ratio for which two-dimen- 
sional values of da/d8 can be used with accuracy has been 
determined for use with the subject report from principles of 
lifting-surface theory and from experiment. An analysis 
based on lifting-surface theory indicates that the da/dS of a 
wing does not vary appreciably with aspect ratio until aspect 
ratio becomes less than approximately 2; for -4<)2, dafd5 
increases rapidly with decreasing aspect ratio until at- -4=0, 
da/d5=l. The discussion of reference 2 concludes that- only 
a lifting-surface theory that accounts for changes in chord 
loading due to flap deflection will give accurate loading pre- 
dictions for low-aspect-ratio wings. 

Low-speed experimental values of dafdS, from, a 7- by 
10-foot wind tunnel (reference 5), are plotted in figure 8 as a 
function of aspect ratio for an unswept and a swept-back 
wing.- Figure 8 indicates the order of variation of da/dS with 
aspect ratio. 

The compressibility considerations iu the present report 
are the same as thpse of reference 1, 2, or 3; namely, that 
compressibility can be taken into account by substituting 
for the geometric values of A, A, and b, the parameters p A, 
Af, and pb, respectively. 

From two-dimensional theory it can be shown that the 
parameter da/dS is independent of Mach number. In the 
present report da/d8 is also considered independent of Mach 
number and, further, is considered independent of aspect 
ratio for -4j> 2, as discussed earlier. However, from the 
linearized compressible-flow theory the aspect ratio becomes 
an effective aspect ratio given by pA which approaches zero 
as M— >1.0. Thus the use of two-dimensional da/dS is re- 
stricted to the Mach number range for which £-4j>2. 

Figure 9 shows the flap-effectiveness comparison between 
theory and experiment for a subsonic range of speed for two 
wings. The theoretical values were obtained using da/dS 
from figure 3 and assuming it independent of Mach number. 
The curves are terminated at the Mach number where 
£-4 = 2. The parallelism of the theoretical and experimental 
curves of figure 9 indicates that dajdb can be used with the 
assumption that it remains independent of Mach number to 
values of pA of the order of two. 

In view of the foregoing discussion, it is believed that 
two-dimensional values of dajdd can be considered independ- 
ent of M as long as £-.4 2:2 in the application of the subject 
procedure., 

COMPARISON OF THEORY WITH EXPERIMENT 

Experimental values of flap effectiveness are compared 
with theoretically predicted values in a correlation diagram 
given by figure 10. Included are the results for a wide range 
of plan forms. Sketches of the plan forms and flaps are 
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drawn about the points of correlation. The theory makes 
use of the curve of figure 3 giving dajdh for a sealed-gap flap 
over a range of deflection of ± 10°. As seen in figure 10, the 
theory is in fair agreement with experiment. 

In general, the theoretical values of figure 10 are too low 
and center approximately about, the minus 10-percent line 
of the figure. A similar effect is shown in figures 8 and 9. 
This almost constant difference between the theoretical and 
experimental results is possibly due to the da/dS values that 
were used in obtaining the theoretical results and which were 
obtained from the —10° to 10° sealed-gap curve of figure 3. 

The theoretical flap effectiveness of the low-aspect -ratio 
wings having partial-wing-chord flaps, given in figure 10, fall 
below the minus 10-percent line. Part of this discrepancy 
is presumably due to using two-dimensional values of da/dS 
since these wings approached the lower limit of aspect ratio 
for which the da/dd concept, has been judged acceptable. 

Comparison of experimental results with theory for a Mach 
number range is shown in figure 9. This figure show's that 
compressibility effects are accurately predicted up to Macb 
numbers corresponding to values of /L4=2. 

CONCLUDING REMARKS 

The determination of symmetric loading due to flap deflec- 
tion for arbitrary wings is shown to be._obtained by the solu- 
tion of four simultaneous equations. The coefficients of the 
simultaneous equations are presented in charts as functions 
of parameters that include wing geometry, compressibility, 
and section lift-curve slope. Integration formulas are given 
for the loading distribution. With these the flap effective- 
ness can be found. 

For the special case of straight-tapered wings, the loading 
distributions and values of flap effectiveness are given in 
chart form for a range of wing plan forms. These charts 
indicate that the effects of plan form on the loading distri- 
bution and lift due to flap deflection are similar to that on 
additional loading. 

On the basis of experimental comparisons, it is concluded 
that the method of the subject report can predict the flap 
effectiveness of wings at subsonic speeds with good accuracy. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Sept. 22, 1950. 


APPENDIX A 

DEVELOPMENT OF EQUATIONS FOR DETERMINATION OF 
SYMMETRIC SPAN LOADING AND LIFT DUE TO FLAP 
DEFLECTION 

ALTERED BOUNDARY CONDITIONS FOR FLAP DEFLECTION 

In the following development the flap chord is considered 
equal in length to Lire wing chord; thus a deflected flap 
results in a discontinuity in the spamvise distribution of 
angle of attack. Spanwise loading for a discontinuous dis- 
tribution of angle of attack can be found through the simul- 
taneous solution method of equation (1). The number of 
equations, however, must be greater than four to obtain 


accurate results and then the computation efforts become 
excessive. (Keference 1 gives the method for an arbitrary 
number of equations.) However, an alternate procedure, 
originally developed in reference 2, provides sufficient accu- 
racy with a reasonable number of equations and this proce- 
dure is the basis of the following development. 

The zero-aspect-ratio theory of reference 3 show's that for 
a given sw'eep angle, as aspect ratio approaches zero, span- 
wise loading becomes independent of plan form. Thus, the 
altered boundary conditions (see Development of Method) 
required for the present theory to give spanwise loading for 
given flap spans can be determined from the loading given 
by reference 3, which applies to all plan forms. These 
altered boundary conditions are obtained from equation (1) 
using the a,» values for zero aspect ratio and the loading dis- 
tribution for zero-aspect-ratio wings given by reference 3 
for given flap spans. 

The zero-aspect-ratio values of a sn are given by figure 1 
for H ,— 0 or are given by the B vn coefficients of reference 1 . 
These a vn coefficients are tabulated as follows: 


TABLE A1 

o,n (FOR m—7 AND .4=0} 


\ v 
n 

1 

2 

3 

4 

, 

10. 4524 

-2. 0720 

0 

-0. 2242 

2 

-3.S2S4 

5.0568 

-L62S4 

0 

3 

0 

-2.3SSS 

4. 3293 

-3. 1543 

4 

2928 

0 

-1.7072 

4. 0000 


With these a r „ coefficients and equation (1), symmetric 
spanwise loading can be found for zero-aspect-ratio wings. 
As a comment on the accuracy of the simplified lifting- 
surface theory for determining additional loading with 
m = 7, the solution of equation (1) with A — 0, a, n values gave, 
to the accuracy of tliree decimal places, the elliptic loading 
distribution characteristic of zero-aspect-ratio wings. 

The flap end is arbitrarily chosen for the present theory 
as the mean value of the spanwise trigonometric coordinates 
of the points at which the boundary condition is applied and 
which bracket the discontinuity or flap spanwise end. Thus 
for m= 7, three flap spans can be defined for both inboard 
and outboard flaps. Let jj/be the flap span, and 8 the span- 
wise point at the end of the flap, then 

r}f = cos 8 for inboard flaps 

ri /= 1 —cos 8 for outboard flaps 

These flap spans are given in the following table: 


TABLE A2 

DEFINED. FLAP SPANS 



Inboard 

Outboard 


' Casa 

I ’ 

II 

hi 

IV 

V 

VI 


$ 

7t 

Sir 

3?r 

3jt 

Sr 

7 r 



16 

16 

16 

16 

16 

16 

_ 


' r 

0. 1951 

0. 5556 

0.8315 

0. 1685 

0.4444 

0. 3049 
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For the flap spans listed in table A2, the exact span load- 
ing distribution for zero-aspect ratio can be found from 
reference 3. With the a„ values listed in table Al and the 
exact values of G lt G 2 , G t , and (? 4 determined from reference 
3, equation (1) gives the boundary condition or twist re- 
quired for the present method to produce the loading dis- 
tribution for each case listed in table A2, that is. 


^•=10.4524 ~ — 3.8284 ^- 2 -0.2928 ^ 

Si Si Si Si 

^=-2.0720 fl+5.6568 ^—2.3888 % 

0 l §1 Ol Oi 

^=-1.8284 -^+4.3296 f^-1.7022 ^ 

01 0] 0 1 Oi 


(Al) 


TABLE A4 

TWIST DISTRIBUTION, ^ 

01 



Inboard 

Outboard 

Case 

I 

n 

m 

IV 

V 

VI 


7 TT 

5x 

3w 

3x 

5x 

7x 


16 

16 

16 

16 

16 

16 


0.1951 

0.5556 

aS315 

0.16S5 

0.4444 

0. S049- 

a\ 

Si 

-.0170 

.0318 

-.0414 

1.0414 

.9682 

1.0170 

cc* 

V 

.0285 

-.0034 

1.0210 

-.0210 

1.0034 

.9715 

03 

-.0136 

.9041 

.954S 

.0452 

.0059 

L0136 

04 

Si 

.9SS4 

.9763 

1.0404 

-.0404 

.0237 

.0116 



-0.2242 ^-3.1548 -^+4.0000 ~ 

o i 5i o i „ 


The spamvise loading distribution due to flap deflection 
is given by (reference 3) : 

Inboard flaps 


g(y) 

Si 


=sin 

7T 



. Q-\-to 

sin 


28 sin ^+(cos <p— cos 6) In 

L 

+ 

. 8~(f 

sm 2 




(cos ^+cos 8) In 
Outboard flaps 


8-{-<p 
COS 2 


6 — ^ 


cos — ~ 


2 

_ 


f (A2) 



With equation (A2), the spanwise loading G u G~, Gt, and G it 

AT < <1P flC 

at- span stations <p=-^> -x-i and ^ or tj = 0.9239, 0.7071, 

0.3827, and 0 can be tabulated for each of the cases given in 
table A2 as follows : 


TABLE A3 


SPANWISE LOADING FOR h=0 


Gu 

St 

Case I 

n 

nr 

IV 

V 

VT 

Gu 

Si 

0.0505 

a 1524 

0.2586 

0.1241 

0. 2303 

0.3322 

Si 

Si 

.1106 

.3467 

.6432 

.0639 

.3604 

.5965 

Gi 

Si 

.2002 

.6903 

.8714 

.0525 

.2336 

.7177 

Gi 

Si 

.4125 

.7970 

.9618 

.0382 

.2030 

.5875 


The twist distribution required for each case is obtained 
with equation (Al) and table A3, and is summarized as 
follows : 


With the twist distribution given by table A4, equation 
(1) can be used to solve for spanwise loading due to flaps 
for any of the six cases and any aspect ratio. 


LIFT DUE TO FLAPS 

The loading distribution resulting from flap deflection is 
distorted to such an extent that the usual quadrature 
formulas are not sufSciently accurate to integrate the loading 
for m = 7. However, the following development results in a 
method which will give the lift due to flaps as a function 
of the four loading coefflcients: 

Lift coefficient is given by the integral 

C l =aJ G(n)dn=A G{G)sm<pd<p 
If it is assumed the loading series 


then 


<7(D = a„ sin m ip 



1 


(A3) 


If the loading series is expanded for <P=^> y’ yp and 


and Gi, G 2 , G 3 , and are obtained in a series of a^’s, then 
multiplying each series through by the first coefficient of the 
series and summing the (Ts gives 



^0.3827 ^TO.7071 ft+0.9239 G 3 + 



Oi 5 — Ui7+ a 3i — O 33 + • • • (A4) 

The high harmonic coefficients can be put as factors of 
the G s . With equation (A4), the lift coefficient of equation 
(A3) becomes : ~ 

®is — u £7 -{- cz 32 — a 3 


C l =A 

0.7071tt 

4 

0.9239 


[0-Ap( H 
(H 


s 0 


1.5068 G, 

U 15 — £l7 + £31 — g 33 

1.506S G 2 

a is — Ol7+®31 — a 33 


~- a * 3 ) (n-r 


1.5068 Gi 


) 


<?2 + 
(?3+ 


77 /l 1 UlS - t_a 17 _ b a 31 a 33 


7T/1 

4 \2 1 


1.5068 Gi 


G< 


] 


(A5) 
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If h n is defined as the coefficients of G n , then the lift 
coefficient is given by 

Cl— A. (h^GiC -]- hiGi) (A6) 


Reference 3 also gives loading due to flaps in scries form 
from which the ratio of (a !5 — a, 7 +a 3 ,— ass) to G„ can be 
evaluated. It is expected that this ratio will not vary appre- 
ciably with aspect ratio. These high harmonic coefficients 
are small but are not negligible for loading due to flaps. The 
a H values from reference 3 are given by the following equa- 
tions (for odd values of p,) : 

For inboard flaps 


Mi= 1 


«i=— (it— 2f?+sin 28) 

7T 




45 , 


For outboard flaps (for the same (?) 
M i = l 

(®l)outboard— 5, (<X i) inboard 

A*i> 1 

(^b,) outboard — (cy) inboard 


Ox, sin Bcospii 8 — cos 6 sin y. x 6) 


Y (A7) 


With equations (A5) and (A7) and the values of table A3, 
the h„ values for each case can be tabulated : 


TABLE A5 

LOADING INTEGRATION FACTORS, h n 


Case- 

• I 

n 

m 

... i . 

IV' 

V 

VI 

hi 

0.2091 

0.2901 

0.2999 

0.3020 

0.3014 

0.3003 , 

hi 

.5541 

.5544 

.5519 

. 5600 

. 5503 . 

.5556 

ht ■ 

.7248 

.7250 

.7252 

.7339 

.7275 

, .7259 

In 

.3922 

.3921 

.3922 ; 

. 40$) 

.3950 

. 3930 . 


INTERPOLATION FUNCTION FOR SPANWISE LOADING DISTRIBUTION 


The method of this report gives spanwise loading due to 
flap deflection at four span stations. Since the complete 
loading distribution is not sufficiently defined by these values, 
the values of loading at other span stations cannot be inter- 
polated accurately by direct use of the interpolation equation 
and table of reference 1 . 

The direct use of the interpolation table of reference 1 can 
be obtained by converting the present loading due to flap 
deflection to a loading approximating additional loading dis- 
tribution. The conversion factors are obtained from values 
of the zero-aspect-ratio theory and given by the ratio of 
equation (A2) to sin <p K for each of the four given span 
stations. Define 


^_ (G/5 i)a-o 
sin <p 


(A 8) 


then, with, table A3 and equation (A8), values of R n have 
been computed and are tabulated for 'each of the flap span 
cases. 


TABLE AO 
R« 



Inboard 

Outboard ’ 


Case 

I 

n 

m 

IV 

V 

VI 

n 

Vt 

V \ 

A A 

0. 195 

0.550 

0.831 

0. 169 

0. 444 

' 

6.805 

I 

0.924 

0. 1318 

0.3981 

a 6757 

0.3243 

0. 6019 

0.SGS2 

2 

' ,707 

.1564 

.4903 

.9096 

— 

.9041 

.5097 

.8430 

3 

;383 

. 2232 

.7471 

.9431 

.5685 

.2529 

.7768 

4 

0 

.4125 ' 

.7970 

. 961*8 

.0382 

.2030 

.5875 


where, for 7^=1, R n = 1. 

The interpolation table of reference 1 gives 
TABLE A7 

INTERPOLATION FACTORS e„ k 



0. 9808 

0.8315 

0. 5556 

6.1951 

\ 





\ ' k 

1 

3 

5 

7 


2 

2 

2 

2 

\ 




V; ■ 

1 

0. 8657 

0. 4904 

, -0. 0970 

; 0.0229 

2 

-.3753 

.7682 

.5132 

-.074.6. 

3 

.2778 

'-.3524 

‘ .7911 

.4157 

4 

-.1275 

.1503 

2250 

.6407 







The interpolated loading at span station k is determined 
from the summation 



where equations (A2) and (A8) provide values of R k which 
are tabulated as follows : 

•• . TABLE A8 
R k 











\.Case 







k 

Vk 

I 

u 

m 

XV 

V 

VI 

1 

2 

• a 981 

0. 1246 

0. 3826 

0.6090 

0.3911 

0. 6176 

0,8754 

3 






. 


2 

.831 

.1421 

.4304 

.8090 

.1910 

.5695 

.8580 

5 








2 

-556 

.1799 

.6253 

.9314 

.0685 

. 

.3746 

.8201 


/ 







7 

2 

.195 

.3324 

.7863 

.9556 

.0442 

.2136 

.0075 


c 

With the R k values of table A8 and the loading ratios -j- 


determined from equation (A9), the loading G k can be found 
13 5 7 

at span stations flr=-> 2 ’ and 2 ( 7 ? ==0 - 981 > 0-831, 0.556 
and 0.195, respectively). 

APPENDIX B 

AIDS TO THE METHOD 

LINEAR ASYMPTOTE OF a,. 

For large values of H„ the a, n functions become linearly 
proportional to H t (see, e. g., fig. 1) . Since this linear charac- 
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teristic appears at relatively low values of H, the simple 
linear relation, between a, n and H, is quite useful. 

The L function of appendix A of reference 1 is multiplied 
by bfc, and the product is linearized with respect to bjc,. 
The result gives 


fb \ T * 2 fb\ , , 1 l\ . . 

\e,) COSA\C,) ' — i?| i7/ bm * 


" „ 2 — 2 • 

V —V 


l . A “V 
-tan ? 


/l+fel tan.O’-l 
if Kn+v / 

2 7j tan A 


for !j>i 


1 1\ . . 2i? . 1 A 

sin A ? — =t+o~ tan -- 

77 — 77 1 nJ V'—V 2 ?j 2 


J-L_ 

\h — ij| 

V i+ ( 


2 v tan A 

A ) (*')_!! 5A+ 

vcos A / Vc,/ t? 


1,1, A 
o — rw-tan^ • • • 
J 7} Ztj Z 


for 7 ?<t 7 


for 77 = 1 


*^CB1) 


= — 2 (— ^— r+tan.0)( -^)+= (1+sinA) . . . 

\cos A ./ \cj t; 

for t? = 0, 

With the values of equation (Bl) substituted into equation 
(A37) of reference 1, the values of a, n for arbitrary sweep 
angle and for high values of H, are obtained. For rn—7, the 
following equation (B2) gives values for a,, as . 


/ Mila-!. 452 tan a'W-K.' 9 15 — 1.751 sin A- 
V cos A / 


0.00.3 tan £+0.088 ( A +0.0016 .»a- A-l ~ 
2 V tan A 


0.052 


0.014 ( - 


V 1+0.0177 tan 2 A— 1~ 

tan A 

VI + 0.1716 tan 2 A— 1 ' 


/(L892 3 _ 793 taEL A H 3.100 + 1.610 sin A- 
\cos A / 

0.008 (an £-0.048 (A±M|^1AA)^ 
0.034 ( Vl+0-0177 toa-A-l \ 

\ tan A / 


0.034 ^1 + 0-01^ tan^A-1 - 


0.052 


Vl +0.1716 tan 2 A-l 


=( °' 2 ^ G +2.476tanA~ N ) 7+-0.0 14-1.727 sin A 

V cos A / 


0.030^1+0.037(^00^=1)+ 


0.125 


Vl +0.0177 tan 2 A— l' 


Vl +0.1716 tan 2 A— T 
tan A 


a 14 =( — — — r- — 1.100 tan A ) Hi— 0.043 + 0.637 sin A — 
V cos A / 


0.060tant-0.0176 ^+Mp^g^ 

2 \ tan A 


3.048^ 


1 + 0.0 177 tan 2 A— 1 


0.08S 


\l +0.1716tan 2 A— 1 


( / 0^31 0 _ 408 tanA N )fJ 2 _ 2. 140 - 0.223 sin A— 

\cos A / 

0.004tan^— 0.063 ^ * t an2A ~ 

2 V tan A I 


0.116 


Vl + 0.0294 tan 2 A-l 


Vl +0.0886 tan 2 1 


0018 /vl+M|Mta^A 

\ tan A / 

!2 =(^i2i + 0.267 tan a) i+,+4.640-0.491sinA+ 

0.01 1 tan^+0.163 ‘ ‘ 

2 \ taDA / 


0.063 


Vl + 0.0294tan 2 A— 1 


0.068 AiL±M88 gtan 2 A-A 
V tan A / 

- a ., s =/^?h£_ 0.986 tan &) Hi-2.226 + 0.821 sinA- 
\cos A / 

nA<0 , A A , „_/Vl +0.0177 tan 2 A— 1\ , 
0.042 tan ^ 0.107^- ^ )+ 

\ tan A / 


0.048 


0.063 


Vl + 0.0886 tan 2 A— 1 


a ,, =( — + 0 . 3 7 7 tan A ) H 2 + 0 . 188 - 0 .551 sir 
" V cosA / 

„ „„„ , A , n ntn ( Vl + 0.0177 tan 2 A— l' 
0.070 tan 2 + 0.048^“ gn . 

Q. 02 2(Vl+M- f»an»A-l V 

\ tan A / 


0.116 


'Vl +0.0886 tan 2 A-l 


/'_0+4t_ 0 _o68 tan aV++ 0. 197+0.207 sinA— 
V cos A / 

0 . 007 tan^- 0 . 116 ( ^ i r 6 r A - 1 )- 

2 \ tan A / 


,/4+Ojputfti 


0.2134 A +0 - 199 g+ A -f 
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(n i — a 

■^—+0.251 tan A ) II 3 - 1.850-0.132sinA+ 


0.020 tail |+0 


0.082 


( yi+o 


\ 


\ tan A / 

)- 


.0886tan 2 A— Li 


tan A 


®33 — ( 


o nc 6 /l+Q ‘ 1993 tanA ~~ l X ) 

\ tan A ) 

0.543 


cos A 


0.232 tan A )ff 3 +4.272-0.388smA- 


A 


0.078 tan 0.197 


71 + 0.1716 tan 2 A— 1 


tan A 


+ 


0.302 


vi + 0.0886 tan 2 A- 


tanA 


a-3 4= 


0.088^ 1+0 ' I993tan2A_ 

/0.156 


tan A 


1 > 

0 


cos A 


0.145tan|+0.088 


0.116 

0.041 


0.428 tanA^// 3 —1.666 + 0.489smA- 

Vl + 0.1716tan il A— 1 \ 
tan A / 

tan 2 A— 1\ 

A ) 


Vl+o. 08801 


y (B2) 


o 


tani 


Vl + 0. 1993 tan 2 A— 


tan A 


0 


a 4 i = 0.032 — ^-— + tan A^)ff 4 — 0.224 
\cos A / 

a 42 = - 0 . 0 9 7 + tan A^ H t - + 0 . 3 54 ( 1 +sin A) 


a 43 =0.371 
a 44 = 0.696 


( 


1 


Vcos A 

(-1 

\cos A 


+ tan A 


)h 4 - 3. 


155 


tan A )i7 4 +3. 313 — 0.687 sin A 


UNEAR SPANWISE DISTRIBUTION OF («),/(«)., 

Equation (3) can be manipulated into a function of various 
parameters involving distributions of wing chord and section 
lift-curve slope. One procedure is as follows: With the con- 
dition that the product of section lift-curve slope and wing 
chord varies linearly spanwise, then 

(kc)i A1.(1 + k< X/*) [} ~ Vr O ' ~i x )] (B3) 


and equation (3) becomes 


H r =d,(pA K ) 


1 +K t \/K, 

2 [l — 17,(1— k,\(k,)] 


where A , is the aspect ratio based on a wing chord equal 
to (kc). Equation (B3) is reduced to terms of two param- 


eters, A k and (a'Vay)X. Expressions of A, in terms of aspect 
ratio for straight-tapered wings, and the distribution of sec- 
tion lift-curve slope can be found. 

For straight-tapered wings 


c, ( 1 + A) 


and since (vc) is linear 
A k 

then 


2 b 


K r C r [1 +(*,/«,-) X] 


A' 


_A(1 + X) 


A'r+ K ( X 


and equation (B3) becomes 


i + 0q/jQx 


(B4) 


L Or + A) J 2 { 1 — n> [1 — (k</(C,)X] } 

The distribution of k for straight-tapered wings is given by 


(kc)„ 1 — r?v [1 — A/+M 

Kr ~c r ~ Kr - 1 — 17,(1— X) 


(B5) 


Equation (B4) is in terms of the two parameters ^~j - J 

and (fc,/<c r )X. Solutions for spanwise loading in terms of these 
two parameters and Ap are valid for the distributions of sec- 
tion lift-curve slope given by equation (B5). Equation (B5) 
indicates that at X = l, k, is a linear function and at X=0, k, is 
a constant. For values of X between 0 and 1, k, is a curve in 
the region between the linear function and a constant. 
Equation (B4) is represented in graphical form by figure- 2 but 

with the ordinate replaced by the parameter- 

and the abscissa by («,/«,) X. 

In summary, for straight-tapered wings having a linear 
distribution (kc)„ the loading and associated aerodynamic 
characteristics can be presented in terms of the parameters 


0A(l+X)/(/q+K(X) 


Ajs, 


§aa±X) ! 

Kr+A'iX 


and (#c ( /i c r )X. 
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Hi 

Fisdre 1.— Influence coefficients, o„, for symmetric spanwise loading plotted as a function of the wing-chord-distribution parameter, lh, for values of the sweep parameter, As, degrees. 

(a) t»l,n»l,'di-0.061. 
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•Figure I.— Continued, (c) y=l, n=3, di=0.06I. 
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Figure 1. — Continued, (i) y=3, n=l, efe= 8.381, 



Figure 1 — 
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Figure I.— Continued* 
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Figube 3. — Variation of lift-effectiveness parameter with flap-chord ratio, — * Average trafling-edge angle about 10°; 21/ = 0.2. CCurves from reference 4.) 
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(a) “-“0; inboard and outboard flaps for all plan forms with the trailing edge at or behind the upstream line of maximum span (reference 3). 
Figure 4. — Spanwise load distribution, -i per radian, due to symmetric flap deflection for straight-tapered wings. - = 1. 

oi c 
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Figure 4 —Continued. 
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(e) X=l.o. 

Figure 5,~Coneluded. 



, per radian 


THEORETICAL SYMMETRIC SPAN LOADING DUE TO FLAP DEFLECTION FOR WINGS AT SUBSONIC SPEEDS 



f a) Effect of aspect ratio for X=0.5. 
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Figure 7.— Flap effectiveness as a function of flap span. 
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Figure 8 .— 1 Variation of experimentally determined ia/is with aspect ratio for wings having flaps of 30 percent or the wing chord normal to the quarter-chord line. 7 /= 1.0, X=0.5. (a) Unswept 

wing; A bin,, iin.-O". (b) Swept wing. 



Figure 9.— Comparison of theoretical and experimental flap-effeetiveness variation with Mach number. Reynolds number a 2X10® at A/=0.2; 4.8X10* at M=0.6. 
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